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Abstract 

Low energy spectra of isotropic quantum dots are calculated in the regime of 
low electron densities where Coulomb interaction causes strong correlations. 
The earlier developed pocket state method is generalized to allow for con- 
tinuous rotations. Detailed predictions are made for dots of shallow confine- 
ments and small particle numbers, including the occurance of spin blockades 
in transport. 
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Much of our present understanding of small quantum dots, with observable discrete level 
structure , concentrates on the regime of relatively high carrier densities where the inter- 
action and charging energy is comparable to the kinetic (Fermi) energy in magnitude 
Similar to real atoms effective single particle orbitals establish a reasonable approximation 
to the electronic states. The spins follow from Hund's rule which is a perturbative result 
though it accords well with experimental findings in small quantum dots at high particle 
densities |Q. 

At lower densities Coulomb interaction is expected to destroy this single particle picture, 
leaving strongly correlated or even crystallized electrons with collective low energy excita- 
tions. While in the homogeneous two-dimensional case rg should exceed Vc = 37 to reach 
this regime (?"s = (^rns)^^/^ measures the ratio between Coulomb and kinetic energy and 
is regulated by the two dimensional carrier density Ug ) , disorder is predicted to reduce this 
value considerably to Tc = 7.5 [§. An even more pronounced reduction of Vc in comparison 
with the homogeneous value is found for the transition into the 'Wigner regime' in quantum 
dots PJTO[]. Careful quantum Monte Carlo (QMC) studies based on the spin sensitivity of 
the density-density correlation function yielded Tc = 4 for parabolic quantum dots [ITU 



Experimentally, this regime has been addressed using capacitance spectroscopy ||Tl[ which 
only probes ground state energies. Non-linear transport behaviour [|l],^|l2|,|l3l has not yet 



been investigated to detect the interesting correlation effects for the low energy excitations. 

Numerical investigations of the low density regime, emphasizing the spin states of rotat- 
ing three electron Wigner molecules, have been carried out for shallow parabolic dots 



Investigations for larger particle numbers have focussed on dots of low symmetry where cor- 
ners in the confining potential or impurities suppress zero modes to delocalize the charges 
in the Wigner regime by so that 'pocket states' can be introduced , which are well suited 
to describe localized charges. The 'pocket states' served as basis to map the spin sensitive 
low energy physics to the one of lattice models of the Hubbard form that account for 
quantum correlations by hopping between nearest places. Applicability of this archetype for 
correlation phenomena has been demonstrated e.g. in quantum dots of polygonal geometry 

i- 

This mapping to a lattice model cannot be carried out straightforwardly if zero modes 
cause charge delocalization which by symmetry actually happens in most experimental quan- 
tum dots. They are fairly well described by an isotropic and in fact parabolic model Mj^T^ 
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H = y-^ + v (1) 

where 

f = >i!E-? + Eg^;^ . (2) 

Here, the effective mass m* and the dielectric constant k are material parameters, and 
Xj {pj) are electron positions (momenta) in two dimensions. This model does not explicitly 
involve spin (as opposed to real atoms spin-orbit coupling is negligible in quantum dots) 
so that all of its eigenstates are simultaneously eigenstates to the square of the total spin 
S"^ to eigenvalues S{S + 1) . The present work extends the pocket state method (PSM) 
to allow for rotational symmetry and compares with results obtained by QMC studies ||1U[| . 
Being based on a recently developed multilevel blocking algorithm to circumvent the 
infamous Fermion sign problem this QMC allows for high accuracy to resolve reliably even 
the low energy spin structure at particle numbers significantly larger than those treatable 
by diagonalizations. 

At low densities the charge carriers form a finite piece of an electron crystal a 



Wigner molecule (WM), that might, classically |]T9|, be arbitrarily oriented. Superposition 
of all of the azimuthal degeneracies leads to an isotropic charge density distribution, as 
required by the symmetry of (0) For analytical progress it is tempting to separate out 
the normal coordinate related with the overall rotation and with total angular momentum 
quantum numbers i (in strictly harmonic confinements i refers to the relative part of 
the Hamiltonian since the center of mass motion just adds integer multiples of ujq to all 
of the eigenvalues and does not affect the spin of any of the states pT| , p2| ). However, the 
remaining normal coordinates then would in general no longer describe identical quantum 
particles obeying Pauli's principle and Fermi (or Bose) statistics but they would correspond 
to linear combinations of such particles. Within the PSM it is crucial to know the result of 
particle permutations in order to assign eventually the correct total spins S to the eigenstates 
and eigenenergies . 

Therefore, we treat all of the possible particle exchanges on equal footing, including 
discrete overall rotations of the WM if they correspond to particle permutations. It depends 
on the geometry of the WM whether rotations by 27r/p with p > 1 leave electron places 
invariant so that the Pauli principle relates £ with S . Such a relationship is well known, 
for instance from the example of solid hydrogen H2, where the even i are necessarily S = 
singlett states while the odd £ are S = 1 tripletts (in this example the spins refer to the 
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protons), the reason being the equivalence of rotations by 180 degrees with the exchange of 
two identical spin-half Fermions. Other examples are discussed in . 

Validity of the PSM requires that the spin sensitive excitation energies A , to be cal- 



culated by this method, should be smaller than charge (plasmon) excitations In the 
absence of continuous symmetries this condition is easily fulfilled at small densities due 
to the almost exponential decay of A ~ exp — ^Jrl . Plasmon energies decrease only 
according to a power law ~ r~^/^ for Coulomb repulsions. With their faster decay 
1/2/ = (27rm* /q°° dr r^n(r))^^ ~ r^"^ (depending on the radial charge density distribu- 
tion n(r) , / is the moment of inertia) the total angular momentum excitations, however, 
still decay faster than the plasmons so that eventually the low energy levels will follow 



only from electron interchanges among the places defining the WM |25], including overall 



rotations by 27r/p , i.e. by processes permuting identical quantum particles [26 



From classical |]T9| as well as from quantum Monte Carlo studies it is known that 
up to < 8 Wigner molecules in the parabolic quantum dots are very symmetric : the 
electrons form one spatial shell ( < 5 ) so that p = N , or one electron occupies the center 
(i.e. p = N — 1). Here we focus on < 6. The method can be generalized straightforwardly 
to larger and more complicated geometries of the WM. 

The transition amplitudes for all possible particle permutations constitute the entries t 
of the pocket state matrix . In the classically forbidden cases t can be estimated within 
the WKB approximation as discussed in [|T^^. The complete potential (0), including 
the interaction, goes into this estimate. Often the most important entries involve only 
two or three adjacent particles, as in quantum dots of polygonal shapes which then 
determine the hopping terms in the equivalent Hubbard model. This is different for the zero 
modes : there a much larger number of particles can be involved into a certain permutational 
transition, such as a rotation by 2it/p in isotropic quantum dots. Corresponding entries tn 
to the pocket state matrix are not of tunneling type and therefore not exponentially small. 
In those cases tn = —p'^/Sn'^I is fixed by the energy constant 1/2/ for rotational excitations 
(/ follows from n{r) ). 

This way all of the relevant entries to the pocket state matrix can be estimated. Its 
diagonalization yields eventually the complete set of low energy eigenvalues. Advantage can 
be taken from the fact that pocket states constitute a faithful representation of the symmetric 
group Sn so that diagonalization can be carried out analytically for small systems, A^ < 4 , 
otherwise numerical help is required. Only irreducible representations [A^/2 + S , A^/2 — S] 
are compatible with Pauli's principle for spin-half Fermions PB| , P7| . This fixes the spin S 
for each eigenvalue. 
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The entries |t| ~ e ^ and |tR,| ^ r^"^ vary differently with the strength of the Coulomb 
interaction so that the ratio t/t^ is a measure for the interaction strength. We use 

1 

y ■= r- > 

1 + t/tn 

ranging from l/(H-(7r^/4)p), since \2t\ cannot exceed the Fermi energy in the non- interacting 
limit, up to unity at strong interactions, y —>■ 1 . 

Figure |I] shows the low energy spectrum versus y for = 3 . Our description is designed 
for evaluating excitation energies, i.e. the differences between the energies of different spin 
states. As expected for weak interactions {y < 0.5 ), the ground state is unpolarized ||22|. A 
transition into the spin polarized ground state S = 3/2 , not found in earlier diagonalization 
studies, is seen above a certain interaction strength which for Coulomb interactions and 
GaAs parameters can be estimated to happen when uq < 0.5 meV |TD|. This result complies 
with the QMC studies and can also be seen when carefully examining Figure 1 of the study 



14 1 of a large quantum dot. We would like to emphasize, that this spin polarization is 



an exact consequence of correlations and not the result of the mean field approximation or 
a magnetic field. In transport experiments, when contacting quantum dots with electron 



reservoirs, it should show up as a 'spin blockade' p8|, since the ground states of A^ = 2 and 
A^ = 3 in sufficiently large quantum dots differ then in spin by more than AS = 1/2 (by 
which entering or escaping single electrons can change spin) since the N = 2 ground state 
(with time reversal symmetry) is always a singlett ||29|| . 

For A^ = 4 (not shown here) we confirm the Hund's rule result of a S = 1 ground 
state, as obtained already in density functional calculations PJ^,|^. New is its persistence 
up to strong interactions. The lowest singlett level S = approaches this ground level 
~ exp —ujq ^^'^ as ujq decreases. The rotationally first excited state i = 1 consists only of 
triplett 5* = 1 levels while the spin polarized level S = 2 belongs to the doubly excited 
rotational state, i = 2 , together with another singlett S = level. 

For A^ = 5 (Figure 0), on the other hand, the polarized state 5 = 5/2 joins the 
unpolarized ground state S" = 1/2 in the lowest rotational level at strong interactions. 
This low energy high spin state makes negative differential conductances in the non-linear 
transport likely, due to the spin blockade [^. Rotationally excited levels consist of S = 1/2 
as well as of S = 3/2 spin states. 

The sixth electron is predicted |jTO|, also classically [jl9|, to occupy the center of a 5-fold 
ring. This complicates the pocket state analysis since new types of pair exchanges appear 
(exchange with the central electron) and also the triple exchange t^ (cyclic permutations of 
three adjacent electrons, including the central one) turns out as important, in accordance 
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with WKB estimates Indeed, the PSM spectra do not compare with the low energy 
levels obtained from QMC unless ^3 is included with a similar magnitude as the pair 
exchanges. 

This demonstrates how our approach complements most favorably the QMC simulations 
for quantum dots which yields abolute values for the many particle energies to high accuracy, 
contrary to the method based on pocket states. Very reliable estimates for the t-parameters 
can be achieved which otherwise would have to be guessed by less trusty approximative 
means. On the other hand, QMC is incomplete for the low energy levels since only the 
lowest eigenenergies to given z-component can be simulated. 

For = 6 and confining energies ujq ^ O.lSmeV (GaAs) we find, with increasing energy, 
the spin sequence 1-0-3-2-1-0-2-1-2-1-1-0. The spin S = 1 indicates another interaction 
induced change in the ground state spin since from the non-interacting levels point of view 

= 6 corresponds to a 'noble gas' configuration implying an unpolarized ground state 
spin S = 0. This result also has to be contrasted with the conjecture S = 2 following 



from a static antiferromagnetic WM |[T^ of pentagonal symmetry. The rotational ground 
state i = includes all possible spin states S* = 1, 0, 3, 2 , with the fully polarized state, 
S = 3 , being lower in energy than the lowest S = 2 state, in accordance with QMC. This 
again suggests possible occurance of negative differential conductances for the transition to 
N = 5. 

In conclusion, generalizing the pocket state method we have developed a description 
for the low density regime in isotropic such as parabolic quantum dots. Low energy levels, 
including spin quantum numbers were determined for < 6 . Detailed predictions are made 
for spin blockades as they should be detectable in linear and non-linear transport through 
shallow quantum dots of confinement energies below ujq <^OAmeV (GaAs). . 
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FIG. 1. Low energy levels using pocket states versus y for A'^ = 3 in units of iR. 
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FIG. 2. Low energy levels using pocket states versus y for N = 5 in units of t^. 
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